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Abstract 

Background: Iodine 125 (1251) seed irradiation can be used as an important supplementary treatment for 
unresectable advanced gastric cancer. Here, we aim to comprehensively elucidate the biological effects induced by 
1251 seed irradiation in human gastric cancer xenograft model by using global expression and DNA methylation 
analyses. 

Methods: The 48 mice bearing NCI-N87 gastric cancer xenografts were randomly separated into 2 groups: sham 
seeds (O mCi) were implanted into the control group (n = 24); 125 I seeds (0.9 mQ) were implanted into the 
treatment group (n = 24). The mitotic index and apoptotic index were evaluated by quantitative morphometric 
analysis of the expression of proliferating cell nuclear antigen (PCNA) and in situ terminal transferase-mediated 
fluorescein deoxy- UTP nick end labeling (TUNEL), respectively. Global gene expression changes induced by 1251 
seed irradiation were analyzed by using Nimblegen Human gene expression array. DNA methylation profile in the 
tumors from control group was investigated with methylated DNA immunoprecipitation (MeDIP) and Nimblegen 
CpG promoter microarrays. The changes in the methylation status of selected genes were further investigated by 
using MeDIP-PCR. 

Results: 1251 seed irradiation suppresses the growth of gastric cancer xenografts in nude mice. PCNA staining and 
tissue TUNEL assays showed that both inhibition of cell proliferation and induction of apoptosis contribute to the 
1251-induced tumor suppression in nude mouse model. Gene expression profiles revealed that the expression levels 
of several hundred genes, many of which are associated with apoptosis or cell cycle arrest, including BMF, MAPK8, 
BNIP3, RFWD3, CDKN2B and WNT9A, were upregulated following 1251 seed irradiation. Furthermore, the up- 
regulation of some of these genes, such as BNIP3 and WNT9A, was found to be associated with irradiation-induced 
DNA demethylation. 

Conclusions: This study revealed that 1251 seed irradiation could significantly induce the up-regulation of 
apoptosis- and cell cycle-related genes in human gastric cancer xenografts. And some of the up-regulation might 
be attributed to 1251-irradiation induced demethylation in gene promoter regions. Collectively, these findings 
provided evidence for the efficacy of this modality for the treatment of gastric cancer. 



* Correspondence: luokaiyuanyunnaniai63.com 
1 Kunming Medical College, Kunming, Yunnan, PR of China 
department of Vascular Surgery; the Second People's Hospital of Yunnan 
Province, Kunming, Yunnan, PR of China 



© 2012 Ma et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BlOlVlGCl CGntfBl Commons Attribution License (http://creativecomrnons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Ma et al. Journal of Experimental & Clinical Cancer Research 2012, 31:61 
http://www.jeccr.eom/content/31/1 /61 



Page 2 of 10 



Background 

Gastric cancer is one of the most frequent cancers in the 
world, and almost of 50% gastric cancer death occurred 
in China [1-3]. Surgery offers the only realistic chance of 
cure; However, many of the patients present with unre- 
sectable tumors at the time of diagnosis. Even with re- 
section, still more than 50% of patients will relapse and 
eventually die of their disease [4,5]. Therefore, non- 
surgical methods have attracted increasing attention. 

In recent years, 1251 implantation has been widely 
used to treat prostate cancer and other tumor types be- 
cause of its ability to offer high precision, little trauma, 
strong lethality, and fewer complications [6-9]. Most re- 
cently, Wang and colleagues applied 1251 implantation 
to treat advanced gastric cancer and found significant 
improvement in clinical symptoms and life quality of 
patients [10]. 

Although the 1251 seed implantation have been suc- 
cessfully applied in clinic, its radiobiological effect and 
underlying molecular mechanism are far from fully 
understood. Recently, Zhuang and colleagues indicated 
that continuous low dose rate irradiation influenced the 
proliferation of cells via MAPK signal transduction. And 
apoptosis was the main mechanism of cell-killing effects 
under low dose rate 1251 irradiation in CL187 cells [6]. 
Besides, Ma and colleagues demonstrated that 1251 ir- 
radiation significantly induced cell apoptosis and inhib- 
ited DNMT1 and DNMT3b expression at 4 Gy in 
pancreatic cancer cells. Thus, the irradiation-induced 
apoptosis and DNA hypomethylation might be two key 
mechanisms underlying the therapeutic effect of low en- 
ergy 1251 seed implantation [11]. However, to date, the 
global molecular changes induced by 1251 irradiation 
have not yet been fully understood. In present study, we 
profiled gene expression in human gastric cancer xeno- 
grafts with microarrays to gain a comprehensive over- 
view of changes induced by 1251 seed irradiation. 

Methods 

Animal model 

The human NCI-N87 cells (3x10 /mouse) were sub- 
cutaneously injected into right dorsal flank of each 
BALB/c-nu/nu nude mouse. After 1-2 weeks of im- 
plantation with tumor cells, when tumors reached ~20- 
30 mm 3, the animals were randomized into control and 
treatment groups (24 animals per group). The 1251 seeds 
(0.9 mCi) were injected into mice in treatment group 
through 18-gauge needles, while ghost seed were 
injected into the mice in control group.The tumor size 
was measured using calipers and the tumor volume was 
estimated by the formula: tumor volume (mm3) = (L x 
W 2) x 1/2, where L is the length and W is the width of 
the tumor. Tumor volumes and body weights were mon- 
itored every 3 days over the course of treatment. The 



tumor weight was measured when the mouse was sacri- 
ficed. Mice were sacrificed after 28 days of treatments 
and tumors were removed and fixed in 10% neutral buf- 
fered formalin for histologic and immunohistochemical 
analyses. All animal procedures were carried out with 
the approval of the Animal Ethics Committee of Kunm- 
ing Medical College. 

Histological analysis of tumors 

Tumors were embedded in paraffin, sectioned at 5 um, 
and stained with H&E (Sigma Aldrich, St. Louis, Mis- 
souri, USA). The mitotic index and apoptotic index 
were assessed by quantitative morphometric analysis of 
proliferating cell nuclear antigen (PCNA) expression 
and in situ terminal transferase-mediated fluorescein 
deoxy-UTP nick end labeling (TUNEL), two established 
markers of proliferation and apoptosis. For PCNA 
localization, formalin-fixed, paraffin embedded sections 
were incubated for 30 min with a mouse monoclonal 
anti-PCNA (Nova Castra Laboratories, Newcastle Upon 
Tyne, UK) at a 1:100 dilution. A peroxidase -conjugated 
antibody to mouse IgG (Abeam Inc., Cambridge, MA, 
USA) was applied followed by diaminobenzidine (Sigma 
Aldrich,St. Louis, Missouri, USA) to localize PCNA in 
the sections. DNA fragmentation was assessed by 
TUNEL, using the Apoptag Peroxidase In situ Apop- 
tosis Detection Kit (Serologicals Corp., Norcross, Ga, 
USA). PCNA- or TUNEL-positive cells were quantified 
in 40 randomly selected high-power fields (x 200) of 
each tissue section. 

RNA extraction 

Total RNA was retracted from tumors using Trizol re- 
agent (Life Technologies Inc., Gaithersburg, Maryland, 
USA) according to manufacturer's instructions. Total 
RNA from each sample was quantified by the NanoDrop 
ND-1000 (NanoDrop Technologies, Montchanin, DE, 
USA) and RNA integrity was assessed by standard de- 
naturing agarose gel electrophoresis. Total RNA from 
one tumor from each mouse (6 tumors per group) was 
used for qRT-PCR analysis, whereas total RNA from 
tumors from four mice per group (12 tumors per group) 
was pooled for each microarray hybridization. 

Microarray analysis 

Microarray analysis of whole-genome gene expression 
profiling was performed using Human 12 x 135 K Gene 
Expression Array (Roche Applied Science, Indianapolis, 
IL, USA). Double-strand cDNA (ds-cDNA) was synthe- 
sized from 5 ug of total RNA using a Superscript 
ds-cDNA synthesis kit (Life Technologies Inc., Gaithers- 
burg, Maryland, USA) in the presence of 100 pmol oligo 
dT primers. ds-cDNA was cleaned and labeled in ac- 
cordance with the NimbleGen Gene Expression Analysis 
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protocol (Roche Applied Science, Indianapolis, IL, USA). 
Microarrays were then hybridized with Cy3 labeled ds- 
cDNA in a hybridization chamber (Roche Applied Sci- 
ence, Indianapolis, IL, USA). After hybridization and 
washing, the slides were scanned using the Axon Gene- 
Pix 4000B microarray scanner (Axon Instruments, 
Union City, CA, USA). Then, the data files were 
imported into Agilent GeneSpring Software (Agilent 
Technologies, Santa Clara, CA, USA) for analysis. Nim- 
bleScan software's implementation of robust multichip 
average offers quantile normalization and background 
correction. The six gene summary files were imported 
into Agilent GeneSpring Software for further analysis. 
Genes that have values greater than or equal to lower 
cutoff of 50.0 in all samples were chosen for data ana- 
lysis. The microarray experiment was independently 
repeated in triplicate for each sample group. Differen- 
tially expressed genes were identified through Fold- 
change and T-test screening. GO analysis and Pathway 
analysis were performed using the standard enrichment 
computation method. 

Real-time polymerase chain reaction (PCR) 

DNase-treated total RNA extracted from each tumor 
sample was reverse transcribed using the Transcriptor 
1st Strand cDNA Synthesis Kit (Roche Diagnostics 
GmbH, Mannheim, Germany). Real-time PCR was per- 
formed for quantitative analysis using SYBR green dye 
(TaKaRa, Tokyo, Japan) on the ABI-Prism 7900HT sys- 
tem (Applied Biosystems, Foster City, CA, USA) accord- 
ing to the protocols recommended by the manufacturer. 
Cycling parameters: pre-denaturation 1 min, 95°C; de- 
naturation 15 s, 95°C; annealing 15 s, 60 °C; extension 
45 s, 72°C, 40 cycles; final extension 5 min, 70°C. The 
fold change was calculated using the 2 -AACt method, 
presented as the fold-expression change in irradiated 
tumors relative to control tumors after normalization to 
the endogenous control, GAPDH. All experiments were 
carried out in triplicate technically. All primers are listed 
in Additional file 1: Table SI. 

Methyl-DNA immunoprecipitation and microarray 
hybridization 

Genomic DNA from tumors from six mice in the con- 
trol group was pooled for Methyl-DNA immunopreci- 
pitation (MeDIP) experiment. MeDIP was performed 
as described previously [12]. Briefly, Genomic DNA 
was sonicated to produce random fragments in size of 
200-600 bp. Four micrograms of fragmented DNA was 
used for a standard MeDIP assay as described. After 
denaturation at 95°C for 10 min, immunoprecipitation 
was performed using 10 |ig monoclonal antibody 
against 5-methylcytidine in a final volume of 500 uL IP 
buffer (10 mmol/L sodium phosphate, pH 7.0), 



140 mmol/L NaCl, 0.05% Triton X-100) at 4°C for 2 h. 
Immunoprecipitated complexes were collected with 
Dynabeads Protein A and M-280 sheep anti-mouse IgG 
(Roche Diagnostics GmbH, Mannheim, Germany) at 
4°C for 12 h, washed with 1 x IP buffer for 4 times, 
treated with Proteinase K at 50°C for 4 h, and purified 
by phenol-chloroform extraction and isopropanol pre- 
cipitation. Immunoprecipitated methylated DNA was 
labeled with Cy5 fluorophere and the input genomic 
DNA was labeled with Cy3 fluorophere. Labeled DNA 
from the enriched and the input pools was combined 
(1-2 ug) and hybridized to a NimbleGen HG18 CpG 
promoter Array (Roche Diagnostics GmbH, Mannheim, 
Germany), which contained all well-characterized 
RefSeq promoter regions [from -800 bp to +200 bp 
transcription start sites (TSSs)]. Array was then washed 
and scanned with Axon GenePix 4000B microarray 
scanner. After normalization, raw data was input into 
SignalMap software (Roche Diagnostics GmbH, Mann- 
heim, Germany) to observe and evaluate the methyla- 
tion peaks. A customized peak-finding algorithm 
provided by NimbleGen was applied to analyze methy- 
lation data from MeDIP-microarray as previously 
described. The algorithm was used to perform the 
modified Kolmogorov-Smirnov test on several adjacent 
probes using sliding windows to predict enriched 
regions across the array. 

MeDIP-quantitative PCR assay 

A MeDIP assay, combined with qPCR, was used to 
evaluate quantitatively the methylation status of candi- 
date genes in the tumors derived from the control and 
1251 treatment groups. MeDIP was performed as 
described above. Purified DNA from the immunopreci- 
pitated DNA complexes and from input DNA was ana- 
lyzed by qRT-PCR on an Applied Biosystems 7900 
Real- Time PCR System. The experiment was per- 
formed in triplicate. The relative changes in the extent 
of gene methylation were determined by measuring the 
amount of detected genes in immunoprecipitated DNA 
after normalization to the input DNA. The primer 
sequences are listed in Additional file 1: Table SI. 

Statistical analysis 

The results of the animal experiments and real-time 
PCR were analyzed using SPSS 13.0 software. (SPSS Inc., 
Chicago, IL, USA) All data were plotted as mean ± 
standard deviation. Student's i-test was used to compare 
values between two independent groups. Differences 
were considered to be significance when p < 0.05. 

Results 

Inhibitory effect of 1125 seed irradiation on the growth 
of gastric cancer 
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The effectiveness of 1251 seed irradiation to inhibit 
the growth of implanted NCI-N87 tumors was exam- 
ined in nude mouse model. There were no significant 
changes in the tumor volumes for the first 10 days of 
the 1251 seed treatment. However, after 13 days, the 
1251-irradiated tumors were much smaller, and sig- 
nificant differences in tumor volumes were observed 
over time between the control and 1251 treatment 
groups (Figure 1A). At day 28, the mice were sacri- 
ficed and tumor weights were measured. Statistical 
difference in the tumor weight was observed between 
the control and treatment groups (Figure IB). All 
these data clearly indicated that 1251 seed implant- 
ation could effectively inhibit tumor growth. Besides, 
the body weights of mice were not affected by the 
1251 irradiation and no obvious radiation-induced 
damage was observed in vital organs of mice (data 
not shown), indicating the safety of 1251 seed 
treatment. 



Effect of 1251 seed irradiation on tumor morphology of 
gastric cancer 

To investigate the effect of 1251 irradiation on the hist- 
ology of NCI-N87 xenografts, tumor sections taken 
from mice in the control and 1251 treatment groups 
were stained with H&E. As shown in Figure 2, the 
histologic appearance of tumors in the control group 
was quite different from that in the 1251 treatment 
group. In the control group, the cancer cells were 
densely arranged with large darkly-stained nuclei and 
obvious karyokinesis. In the treatment group, large nec- 
rotic regions were observed around the 1251 seed. The 
cancer cells adjacent to the necrotic region were loosely 
arranged with condensed nuclei and reduced eosino- 
philic cytoplasm. These results indicated that 1251 seed 
implantation caused growth inhibition of cancer cells in 
NCI-N87 xenografts. 



Effect of 1251 seed irradiation on cell apoptosis and 
mitosis of gastric cancer 

To quantitatively compare the mitotic and apoptotic 
index of tumors treated with 1251 seed irradiation, 
immunostainings for PCNA and TUNEL assays were 
performed. As shown in Figure 3A, the number of 
PCNA- positive cells in the 1251 treatment group was 
obviously less than that of control group. And the mi- 
totic index was significantly decreased in irradiated 
tumors as compared to the tumors in the control group 
(Figure 3B). In contrast to the mitotic index, 1251- 
irradiated tumors showed increased numbers of apop- 
totic cells with condensed and irregularly shaped nuclei, 
staining positively for TUNEL (Figure 3 C). the apoptotic 
index was significantly increased in the 1251 treatment 
group as compared to the control group (Figure 3D). 

Identification of genes induced by 1251 seed irradiation 

Gene expression microarrays were used to characterize 
the gene expression changes in NCI-N87 tumors be- 
tween the 1251 treatment group and control group. 
When the Fold Change (FC) is set > 1.3 and the 
p value at < 0. 05, we found that 544 genes were 
induced by 1251 seed irradiation, while 368 genes were 
repressed (Additional file 2: Table S2). To identify the 
biological processes that were induced by 1251 seed ir- 
radiation, Gene Ontology (GO) functional analysis was 
performed. GO terms for biological processes were 
assigned to these differential genes and this procedure 
was essential to provide an overview of the effect of 
1251 seed implantation in NCI-N87 xenografts. Accord- 
ing to GO functional analysis, the categories cell cycle, 
induction of apoptosis, cell division and growth were 
most significantly overrepresented among the 125- 
irradiation induced genes (Additional file 3: Table S3). 
And many of these genes are critical pro-apoptotic 
molecules or genes associated with cell cycle arrest, 
such as MAPK8, BNIP3 and CDKN2B (Table 1). Then, 
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Figure 1 Effect of 1251 seed irradiation on the tumor volume and tumor weight. (A) Tumor volumes. (B) Tumor weight. Data are the 
mean ± SD and analyzed by the Mann-Whitney U test (*: P < 0.05). 
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we employed DAVID software on the basis of the 
KEGG pathway map to further investigate key path- 
ways linked to these genes. Our analysis yielded 11 
pathways, including cell cycle pathway and several 
pathways associated apoptosis and cell cycle arrest, 
such as MAPK and TGF-beta signaling pathways 
(Additional file 4: Table S4). 

To independently confirm the microarray results, real- 
time RT PCR was performed on samples from BALB/c 



mice that had been exposed to the same experimental 
conditions that were used in microarray assay. The rela- 
tive expression levels of six genes — BMF, MAPK8, 
BNIP3, RFWD3, CDKN2B and WNT9A— were assayed 
in irradiated and non-irradiated tumors. There was a 
close correlation between microarray data and qRT-PCR 
data (Figure 4), indicating the accuracy of our micro- 
array data and the significant induction in the expression 
of selected genes following irradiation. 
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Figure 3 PCNA and TUNEL analyses for tumor tissue. (A) Tumor sections immunostained with an antibody against PCNA revealed that there 
were many strongly positive nuclei in control tumor tissues, whereas such nuclei were rare in tumor tissues of 1251 treatment group. (B) 
Quantification of PCNA staining showed mitotic index of 1251-implanted tumor was much lower than that of control group (*: P < 0.05). (C) 
Apoptosis of tumor tissues in different groups were evaluated by TUNEL assays, which showed that 1 25I treatment induced a significant 
enhancement of apoptotic cells in contrast to control group. (D) Quantification of TUNEL assay showed that apoptosis index of 1251-implanted 
tumor was much higher than that of control group (*: P < 0.05). 
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Table 1 1251-irradiation induced genes associated with apoptosis and cell cycle arrest 



GENE NAME 



DESCRIPTION 



Fold change 



P value 



FDR 



Pro-apoptotic genes 

BNIP3 

MAPK8 

BCL2L1 1 

AKT1 

BMF 

P2RX7 

TNFRSF10B 

APH1A 

TRAIP 

JAK2 

TRIM35 

ITSN1 

TAP2 

ACVR1B 

Genes associated with 

CDKN2B 

RFWD3 

HUS1 

PMP22 

CDC25C 

WNT9A 

CCNG2 

CDC7 

C12orf32 



BCL2/adenovirus E1B 19 kDa interacting protein 3 

mitogen-activated protein kinase 8 

BCL2-like 1 1 (apoptosis facilitator) 

v-akt murine thymoma viral oncogene homolog 1 

Bcl2 modifying factor 

purinergic receptor P2X, ligand-gated ion channel, 7 

tumor necrosis factor receptor superfamily, member 10b 

anterior pharynx defective 1 homolog A (C. elegans) 

TRAF interacting protein 

Janus kinase 2 (a protein tyrosine kinase) 

tripartite motif-containing 35 

intersectin 1 (SH3 domain protein) 

transporter 2, ATP-binding cassette, sub-family B (MDR/TAP) 
activin A receptor, type IB 
cell cycle arrest 

cyclin-dependent kinase inhibitor 2B (pi 5, inhibits CDK4) 
ring finger and WD repeat domain 3 
HUS1 checkpoint homolog (S. pombe) 
peripheral myelin protein 22 
cell division cycle 25 C 

wingless-type MIVTTV integration site family, member 9A 
cyclin G2 

CDC7 cell division cycle 7 (S. cerevisiae) 
chromosome 12 open reading frame 32 



2.1 
1.7 
1.9 
1.4 
1.5 
1.4 
1.4 
1.4 
1.4 
1.6 
1.3 
1.5 
13 
1.6 

1.3 
1.3 
1.4 
1.5 
1.5 
1.6 
1.5 
1.4 
1.5 



0.045 

0.017 

5.39E-04 

0.028 

0.005 

0.004 

0.003 

0.010 

0.032 

0.011 

0.018 

0.020 

0.024 

0.009 

0.034 
0.040 
0.017 
0.042 
0.017 
0.048 
0.003 
0.016 
0.002 



0.050 
0.047 
0.036 
0.049 
0.040 
0.040 
0.038 
0.039 
0.046 
0.045 
0.046 
0.046 
0.048 
0.046 

0.049 
0.050 
0.047 
0.050 
0.047 
0.050 
0.038 
0.049 
0.033 



Collectively, these data indicated that many critical 
molecules and pathways associated with apoptosis and 
cell cycle arrest were activated by 1251 seed irradiation 
in NCI-N87 xenografts, thereby highlighting their im- 
portant roles in 1251 irradiation-induced inhibition of 
tumor growth. 

DNA methylation analysis of 1251-irradiation induced 
genes 

Aberrant DNA hypermethylation is commonly asso- 
ciated with cancer. The Dnmtl DNA methyltransferase 
is responsible for maintenance of the DNA methylation 
pattern. Consistent with previous study [11], significant 
decrease of DNMT1 expression was observed in our 
array data, and this result was validated via the real time 
RT-PCR (Figure 5A). These data suggest that DNA 
demethylation might be involved in 1251-induced tumor 
suppression. Because promoter demethylation is asso- 
ciated with gene re-activating, we focused our attention 
on the 1251 irradiation-induced genes by coupling global 



gene expression and methylation profiles. The genes 
with promoter hypermethylation in the non-irradiated 
tumors were indentified with MeDIP-chip analysis (Add- 
itional file 5: Table S5). Among them, we identified 20 
genes whose expression was significantly upregulated in 
the irradiated tumors as compared to the non-irradiated 
tumors (Table 2). Thus, we speculated that the expres- 
sion levels of these 20 genes might be modulated via the 
promoter demethylation induced by 1251 irradiation. 
Notably, several of these genes were associated with 
apoptosis or cell cycle arrest, such as BNIP3, WNT9A 
and GSG2. To confirm our hypothesis, methylation sta- 
tus of these three genes was examined with MeDIP-PCR 
assay in the treatment and control groups. As shown, 
BNIP3 and WNT9A in 1251 treatment group displayed 
lower levels of methylation status compared with control 
group (P < 0.05), which decreased to 50.9% and 41.0%, 
respectively (Figure 5B). Meanwhile, the expression 
levels of BNIP3 and WNT9A were significantly upregu- 
lated in the treatment group (Figure 5C). These data 
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Figure 4 Quantitative RT-PCR validation for differential genes 
in microarrays. (A) Relative mRNA expression of 6 selected genes 
in microarrays (B) Validation of relative mRNA expression of selected 
genes with qRT-PCR. The significance of the varieties between the 
control group (solid bar) and 125! treatment group (hollow bar) was 
analyzed through student's t-f test. (*: P < 0.05). 



indicated that some apoptosis- and cell cycle-related 
genes could be activated by the demethylation of their 
promoters, which were induced by 1251 seed irradiation. 

Discussion 

Several recent studies have suggested that apoptosis and 
cell cycle arrest may have important roles in the thera- 
peutic effects of the continuous low-energy 1251 irradi- 
ation. However, the comprehensive evidences on this 
topic, especially in molecular levels, still lack. In this 
study, microarray analysis of human gastric cancer xeno- 
grafts exposed to 1251 seed irradiation were performed 
to gain insight into the mechanisms underlying the bio- 
logical effects of 1251 irradiation. 

N87 gastric cancer cells were implanted into the nude 
mice to create the xenograft animal model. The growth 
curves of tumors indicated that irradiation induced 



significant tumor growth inhibition. By observing H.E. 
staining slides, large numbers of apoptotic cells were 
observed in gastric cancer receiving 1251 seeds implant- 
ation. Furthermore, the mitotic index and apoptotic 
index were assessed by quantitative morphometric ana- 
lysis of PCNA expression and TUNEL, respectively. In 
our work, a declined mitotic index and increased apop- 
totic index were discerned in 1251 treatment group com- 
pared with control group, which suggests that 1251 seed 
irradiation can restrain tumor growth and lead to apop- 
tosis of cancer cells. 

Next, we use microarray gene expression profile ana- 
lysis to study the mechanism of irradiation-mediated 
prevention of gastric tumors. To our knowledge, this is 
the first investigation to use microarray technology to 
study the role of 1251 seed irradiation in cancer treat- 
ment. At 28 days following 1251 seed irradiation, the 
nude mice were sacrificed and gene expression was pro- 
filed in the xenografts by using gene expression microar- 
rays. We found that the expression levels of 544 genes 
were significantly induced by 1251 seed irradiation. 
Interestingly, among the irradiation-induced genes, 
many are involved in cell cycle, apoptosis and cell div- 
ision. The main pathways linked to these genes were 
further investigated by KEGG analysis and several 
apoptosis- or cell cycle-related pathways, such as MAPK 
and TGF-beta pathways, were clearly indentified. Then, 
the expression of 6 genes (BNIP3, MAPK8, BMF, 
RFWD3, CDKN2B and WNT9A), which were asso- 
ciated with apoptosis or cell cycle arrest, was further 
validated via real time PCR analysis (Figure 3). BNIP3 
(BCL2/adenovirus E1B 19 kDa interacting protein 3) is 
a proapoptotic member of the Bcl-2 family and its mu- 
tation and dysregulation might play a role in gastric car- 
cinoma development [13]. Recent study revealed that 
BNIP3 might play a role in enhancement of radiother- 
apy efficiency, and its expression might have a synergis- 
tic effect on radiation treatments [14]. MAPK8 
(Mitogen-activated protein kinase 8) is a member of the 
MAP kinase and JNK family. This gene is involved in 
UV radiation-induced apoptosis, which is thought to be 
related to the cytochrome c-mediated cell death path- 
way [15]. BMF (Bcl-2-modifying factor) is a Bcl-2 family 
member bearing only the BH3 domain and an essential 
inducer of apoptosis [16]. BMF contributes to enhancing 
effects on apoptosis after ionizing radiation [17]. 
RFWD3(ring finger and WD repeat domain 3) is an E3 
ubiquitin ligase that positively regulates p53 levels and 
regulates Gl Checkpoint in Response to ionizing radi- 
ation [18]. CDKN2B (Cyclin-dependent kinase 4 inhibi- 
tor B) belongs to a family of cyclin-dependent kinase 4 
inhibitors (INK41) and controls cell proliferation during 
the Gl phase of the cell cycle [19]. The expression of 
this gene was found to be dramatically induced by TGF 
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Figure 5 Effects of 1251 irradiation on gene methylation and mRNA expression in xenografts. (A) Relative expression of DNMT1 was 
detected using qRT-PCR. (B) Effects of 1251 irradiation on gene methylation of BNIP3 and WNT9A in xenografts assayed by MeDIP-PCR. BNIP3 and 
WNT9A in treatment group displayed lower level of methylation when compared with control group. (C) Relative expression of BNIP3 and 
WNT9A was detected using qRT-PCR. Data are expressed as the mean ± SD of 6 samples. The significance of the varieties between the control 
group and 1251 treatment group was analyzed through student's t-t test. (★: P < 0.05). 



beta, which suggested its role in the TGF beta induced 
growth inhibition [20]. WNT9A is a member of the 
WNT gene family and over-expression of t human 
Wnt9a induced cell-cycle arrest at Gl/S boundary [21]. 

Consistent with previous study [11], we found signifi- 
cantly decreased expression level of DNMT1 in the 



irradiated xenografts. DNMT1 is responsible for precise 
duplicating and maintaining the pre-existing DNA 
methylation patterns after replication [22]. Therefore, it 
is reasonable to speculate that DNA hypomethylation 
induced by 1251 irradiation might be associated with 
tumor growth inhibition. By coupling data derived from 



Table 2 The irradiation-induced genes with promoter hypermethylation in the non-irradiated tumors 



GENE_NAME 


DESCRIPTION 


Fold change 


Regulation 


P-value 


FDR 


DRD5 


dopamine receptor D5 


1.4 


up 


2.85E-04 


0.03 


PFN2 


profilin 2 


1.4 


up 


0.021 


0.05 


SKI 


v-ski sarcoma viral oncogene homolog (avian) 


1.6 


up 


0.005 


0.04 


WNT9A 


wingless-type MMTV integration site family, member 9A 


1.6 


up 


0.048 


0.05 


CXorf12 


chromosome X open reading frame 12 


2.0 


up 


0.012 


0.05 


BNIP3 


BCL2/adenovirus E1B 19 kDa interacting protein 3 


2.0 


up 


0.045 


0.05 


CHST10 


carbohydrate sulfotransferase 10 


2.2 


up 


0.010 


0.05 


PNMA1 


paraneoplastic antigen MA1 


1.3 


up 


0.001 


0.04 


C18orf55 


chromosome 18 open reading frame 55 


1.4 


up 


0.009 


0.05 


TRAK2 


trafficking protein, kinesin binding 2 


1.3 


up 


0.047 


0.05 


LRRC49 


leucine rich repeat containing 49 


1.5 


up 


0.041 


0.05 


EPB41L4B 


erythrocyte membrane protein band 4.1 like 4B 


1.4 


up 


0.027 


0.05 


USP31 


ubiquitin specific peptidase 31 


1.5 


up 


0.021 


0.05 


GSG2 


germ cell associated 2 (haspin) 


1.6 


up 


0.035 


0.05 


ATAD1 


ATPase family, AAA domain containing 1 


1.3 


up 


0.006 


0.04 


MGC16385 


hypothetical protein MGC16385 


1.4 


up 


0.046 


0.05 


TCEB3C 


transcription elongation factor B polypeptide 3 C (elongin A3) 


2.0 


up 


0.006 


0.04 


LONRF1 


LON peptidase N-terminal domain and ring finger 1 


1.4 


up 


0.014 


0.05 


SAMD1 1 


sterile alpha motif domain containing 1 1 


1.4 


up 


0.031 


0.05 


SLC35E2 


solute carrier family 35, member E2 


1.3 


up 


0.027 


0.05 
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gene expression microarrays with that of MeDIP-chip, 
we found 39 candidate genes whose expression might be 
activated by 1251-induced DNA demethylation. Notably, 
several of the candidates are pro-apoptotic molecules or 
genes associated with cell cycle arrest, such as BNIP3, 
WNT9A and GSG2 (Serine/threonine-protein kinase 
haspin). The promoter demethylation of BNIP3 and 
WNT9A after receiving 1251 irradiation was then suc- 
cessfully validated with MeDIP-PCR. DNA methylation 
of the BNIP3 promoter was mediated by DNMT1 via 
the MEK pathway [23]. Aberrant methylation of BNIP3 
was also detected in 66% of primary colorectal and 49% 
of primary gastric cancers. Epigenetic alteration of 
BNIP3 is a frequent and cancer-specific event, which 
suggests that inactivation of BNIP3 likely plays a key role 
in the progression of some gastrointestinal cancers and 
that it may be a useful molecular target for therapy [24]. 
Methylation of WNT9A promoter occurs frequently in 
primary colon cancers and WNT9A hypermethylation in 
cancer points to its possible role as a tumor suppressor 
gene [25]. 

This study provides first demonstration for the global 
induction of apoptotic and cell cycle-related genes by 
1251 seed irradiation. And some of the induction may be 
mediated by the irradiation-induced DNA demethyla- 
tion, suggesting that 1251 seed irradiation affects genes 
associated with apoptosis and cell cycle arrest in both 
transcriptional and epigenetic levels. Collectively, these 
data provide an explanation for the tumor inhibitory 
effect of 1251 seed implantation and emphasize the im- 
portant roles of apoptosis and cell cycle arrest under- 
lying the efficacy of this modality. 
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